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(U)  LIST  OF  SYMBOLS 


(U)  Geophysical  Parameters 

g  =  Acceleration  of  gravity, 
p  =  Atmospheric  pressure, 
p  =  Density  of  air. 

9  «=  Virtual  temperature,  degrees 
Ke Ivin . 

R  -  Gas  constant  of  dry  air. 

D  =  D- value,  equal  to  altitude  of  a 
given  pressure  in  actual  atmos¬ 
phere  minus  altitude  of  same 
pressure  in  standard  atmosphere. 

z  >»  Altitude  above  mean  sea  level. 

(o  =  Non-dimensional  pressure 

p/2P 

variable  =  J  exp  £  din  £ 

-00 

(u)  Re-Entry  Vehicle  Parameters 

W  =  Weight  of  re-entry  vehicle. 

=  Drag  coefficient. 

=  Hypersonic  drag  coefficient, 

A  =  Cross  section  normal  to  flow. 

P  =  sin  Og  W/CpA  =  atmospheric  pres¬ 
sure  at  peak  deceleration. 


(U)  Trajectory  Parameters 

t  =  Time . 

V  =  Speed  of  re-entry  vehicle. 

Vg  =  Speed  at  300,000  feet. 

a  =  Path  angle,  i.e.,  angle  be¬ 
tween  horizontal  plane  and 
trajectory  tangent. 

Og  =  Path  angle  at  300,000  feet. 

I  =  Time  integral  of  drag  deceler¬ 
ation. 


(U)  Mathematical.  Subscripts, 
Superscripts,  and  Indices 

a(x)  =  Standard  deviation  of  x 
=  [ (x  -  x)8]^ 

x  =  Indicates  average  of  x  over 

sample  representing  locale 
and  month. 

r(x,y)  =  Coefficient  of  linear  corre¬ 
lation  between  x  and  y. 

(  )  =  Denotes  average  with  respect 

to  <0  along  a  vertical 

“>B 

=  (a^  -  cDg)-1  J  (  )  dto 


6 (  )  =  Indicates  deviation  from 

conditions  in  standard  at¬ 
mosphere  and/or  from  condi¬ 
tions  of  ideal  behavior  of 
fuzing  system  components. 

(  )*  =  The  asterisk  indicates  that 

the  quantity  represents 
standard  atmosphere  condi¬ 
tions,  and/or  programmed  be¬ 
havior  of  fuzing  system  com¬ 
ponents  . 

(  )B  =  Refers  to  burst. 

(  )  =  Refers  to  signal. 

E7(x)  =  Exponential  integral 

x 

=  J'  exp  £  din  £ 


(U)  Fuzing  System  Parameters 

T  =  Time  elapsed  between  signal  and 
burst. 

6T  =  Error  of  timer. 

e*  =  Relative  error  of  integrating 
accelerometer. 
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(U)  LEAST  DISPERSION  OF  THE  MARK  12  RE-ENTRY  VEHICLE 


1.  (U)  Introduction 

(U)  The  ever-changing  atmosphere  causes  a  dispersion  of  the  trajectories 
of  ballistic  re-entry  vehicles.  Computer  programs  are  In  existence  which  can 
compute  with  high  precision  the  trajectories  through  a  given  atmosphere.  In 
this  manner  the  effect  of  the  variable  atmosphere  on  the  trajectory  dispersion 
may  be  calculated  by  solving  the  re-entry  equations  for  an  array  of  atmospheres 
which  represent  the  atmospheric  variability. 

(U)  The  Monte  Carlo  method  used  in  conjunction  with  numerical  Integration 
by  electronic  computers  is,  howevei)  expensive .  The  trajectory  program  present¬ 
ly  used  at  BMD/AFSC  costs  about  $100  per  trajectory.  As  an  example,  the  cost 
of  processing  the  meteorological  data  used  in  this  report  by  Monte-Carlo- ing 
the  trajectories  sounding  by  sounding  would  have  amounted  to  1.6  million 
dollars.  The  data  consists  of  four  years  of  radiosondes  between  sea  level  and 
100,000  feet,  two  soundings  per  day  for  two  locations  and  for  one  month  in  the 
summer  and  one  month  in  the  winter,  in  all,  about  1000  soundings.  A  16-point 
change  in  the  trajectory  program,  varying  the  burst  altitude  and  one  fuzing 
parameter  (the  so-called  signal  altitude,  as  defined  later  in  the  report), 
brings  the  costr  up  to  the  quoted  figure. 

(U)  Realizing  that  the  Monte  Carlo  method  is  not  the  best  way  to  make  use 
of  the  vast  collection  of  upper-air  data  already  available  in  the  meteorologi¬ 
cal  archives  in  processed  form,  the  Air  Weather  Service  has  developed  an  alter¬ 
nate  method  of  assessing  the  environmental  effects.  The  method,  which  may  be 
referred  to  as  the  perturbation  method,  has  been  successfully  applied  to  several 
re-entry  problems  [1]  [2]  [3].  Appendices  A  and  B  to  this  report  contain  a 
mathematical  precis  of  the  application  to  the  Mark  12  Re-Entry  Vehicle. 

(U)  The  following  sections  define  the  problem,  present  the  findings,  and 
give  conclusions  and  recommendations  for  the  design  of  the  fuzing  system  of  the 
Mark  12  Re-Entry  Vehicle. 

2.  (u)  The  Fuzing  Problem  of  the  Mark  12  Re-Entry  Vehicle 

(U)  The  object  of  a  fuzing  system  is  to  make  the  warhead  detonate  at  the 
desired  point  on  the  trajectory. 

(S)  The  Mark  12  Re-Entry  Vehicle  uses  an  inertial  fuzing  system.  An  ac¬ 
celerometer  measures  the  drag  deceleration  and  Integrates  it  over  time.  When 
a  predetermined  value  of  the  Integral  is  reached,  a  timing  device  starts  count¬ 
ing  time  at  an  altitude  which  in  this  report  is  referred  to  as  the  signal 
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altitude .  After  a  predetermined  time  interval  has  elapsed,  the  detonater  is 
activated. 

(u)  Since  the  drag  deceleration  depends  on  the  atmospheric  density,  vari¬ 
ability  in  the  latter  will  cause  a  dispersion  of  the  detonation  point  along  the 
trajectory,  with  a  corresponding  dispersion  in  burst  height  and  range. 

(S)  The  errors  of  the  fuzing  components  will  aloo  contribute  to  the  dis¬ 
persion.  These  are  mainly  the  error  in  the  output  of  the  Integrating  acceler¬ 
ometer  and  the  timing  device. 

(U)  The  atmospheric  variability  and  the  errors  of  the  fuzing  components 
combine  to  give  the  burst  dispersion. 

(S)  The  problem  which  this  report  deals  with,  is  to  define  the  signal 
altitude  which  gives  the  least  dispersion  when  the  system  is  applied  in  a  given 
location  and  season  and  has  a  stated  level  of  performance  of  the  fuzing  compo¬ 
nents  . 

(S)  The  analysis  of  the  performance  of  the  Mark  12  Re-Entry  Vehicle  is 
contained  in  Appendices  A  and  B.  The  result  of  this  analysis  is  that  the  burst- 
height  error  dZg  (positive  for  a  too-high  burst)  can  be  written  as  a  linear 
combination  of  four  sources  of  errors  (equation  A-26): 

dzB  =  Dg  +  MS  -  Ne  -  Q6T 

(S)  The  quantities  in  the  above  expression  are: 

a.  (S)  M,  N,  and  Q  are  positive  constants  for  a  given  burst  altitude 
and  signal  altitude.  Their  variation  with  the  altitudes  of  burst  and  signal 
are  illustrated  in  Figures  4,  5,  and  6. 

b.  (S)  Dg  is  the  so-called  D-value  at  the  programmed  burst  altitude. 

.  .  A 

c.  (S)  9  is  the  average  departure  of  temperature  from  nominal  values 
in  the  layer  from  signal  to  burst. 

d.  (S)  e  is  the  relative  error  (fraction  of  total  value)  of  the  inte¬ 
grating  accelerometer. 

e.  (S)  dT  is  the  error  of  the  timing  device. 

(U)  By  means  of  a  statistical  processing  outlined  in  Appendix  B,  this 
expression  has  been  applied  to  two  locations  in  the  U.S.S.R.  for  one  month  in 
summer  and  one  in  winter.  The  results  are  given  in  the  following  section. 

3.  (U)  Burst-Height  Errors  for  Two  Locations  in  the  U.S.S.R.  in  Summer  and  in 

Winter 

(U)  Figure  1  shows  the  contributions  to  the  variance  of  the  burst-height 
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error  at  station  23552  during  January.  This  station  is  TARK-SAEE  (64*55'N- 
77#49'E). 

(U)  The  variance  is  due  to  different  sources.  Panels  A,  B,  D,  and  E  of 
Figure  1  give  the  contributions  of  the  four  sources  in  the  order  they  are 
listed  in  the  previous  section. 

(S)  Panel  C  shows  the  contribution  due  to  the  correlation  of  mean  tempera¬ 
ture  and  D-value.  It  is  assumed  that  the  co- variance  between  the  errors  of  the 
fuzing  components  (accelerometer,  timer)  and  between  these  and  the  meteorologi¬ 
cal  variables  are  insignificant.  There  is  no  reason  to  assume  otherwise. 

(S)  It  is  obvious  from  Figure  1  that  the  timing  error  of  one  millisecond 
has  no  practical  influence.  The  other  contributions  are  of  similar  magnitude, 
though  as  is  obvious,  the  meteorological  contributions  to  the  variance  dominate. 

(S)  Panel  F  shows  the  composite  error.  A  locus  of  minimum  dispersion  de¬ 
fines  the  signal  altitude  which  results  in  the  least  dispersion  for  any  given 
burst  altitude  between  zero  and  20,000  feet.  For  programmed  bursts  between 
zero  and  7000  feet  the  optimum  signal  altitude  is  about  7000  feet  above  the 
burst  altitude.  For  bursts  between  7000  feet  and  12,000  feet  the  optimum 
thickness  of  the  layer  from  signal  to  burst  shrinks  gradually,  reaching  zero 
at  12,000  feet.  This  means  that  the  timer  should  be  dispensed  with  above 
12,000  feet  (if  practical). 

(S)  As  this  example  shows,  the  exact  location  of  the  signal  altitude  in 
relation  to  the  burst  altitude  is  not  very  critical.  The  minimum  dispersion 
(one  sigma)  is  about  430  feet  for  this  location  and  month  for  bursts  between 
zero  and  5000  feet  altitude.  Moving  the  signal  altitude  up  or  down  3000  feet 
from  the  optimum  position  would  increase  the  dispersion  by  only  20  to  30  feet. 

(u)  Figure  2  shows  the  burst-height  dispersion  (one  sigma)  for  the  two 
stations  during  January  and  July:  23552,  TARK-SALE,  U.S.S.R.  (64#55'N, 

77°49'E)  and  38457,  TASHKENT,  U.S.S.R.  (4l°l6'N,  69#l6'E).  The  results  are 
based  on  four  years  of  data,  1955-1958,  or  about  250  soundings  for  each  loca¬ 
tion  and  month. 

(S)  In  the  four  cases  the  loci  of  minimum  dispersion  are  all  similarly 
located.  It  appears  that  latitudinal  and  seasonal  aspects  do  not  greatly  in¬ 
fluence  the  position  of  the  minimum-dispersion  curve. 

(U)  Figure  3  shows  the  minimum-dispersion  curves  superimposed  on  a  grid 
indicating  approximately  the  time  of  flight  between  signal  and  burst.  This 
diagram  permits  one  to  compare  the  minimum-dispersion  curves  with  any  program¬ 
ming  restrictions  placed  on  the  timer. 

(S)  Observe  that  a  timer  setting  of  three  seconds  for  bursts  between  zero 
and  5000  feet  and  a  setting  of  two  seconds  for  bursts  between  5000  feet  and 
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10,000  feet  will  indicate  the  least  dispersion  fairly  well. 

(U)  Figure  3  carries  a  scale  in  the  right  margin  which  shows  approximately 
the  output  of  the  integrating  accelerometer  in  g  seconds.  Any  programming  re¬ 
strictions  due  to  bounds  on  the  output  of  the  accelerometer  may  be  evaluated  by 
means  of  this  scale. 

(S)  These  results,  while  calculated  here  for  four  specific  atmospheres, 
have  a  more  general  validity.  The  D-value  and  the  mean  temperature  vary  in 
such  a  fashion,  irrespective  of  season  and  location,  that  the  signal  altitude 
which  gives  the  least  burst-height  dispersion,  is  found  6000  to  7000  feet  above 
the  burst  altitude  for  bursts  below  10,000  feet. 

(S)  The  results  above  are  predicated  on  a  relative  error  of  the  integrat¬ 
ing  accelerometer  of  0.5  percent.  A  larger  value  will  tend  to  locate  the 
least-dispersion  curve  higher  up. 

4.  (u)  Conclusions  and  Recommendations 

(U)  As  a  result  of  this  study,  the  following  conclusions  have  been  reached 
and  recommendations  are  made: 

a.  (S)  The  Mark  12  Re-Entry  Vehicle  will  have  errors  in  the  burst 
altitude  (and  range)  which  are  due  to  the  errors  of  the  fuzing  components  and 
to  the  departure  of  the  atmospheric  structure  from  the  particular  atmosphere 
used  in  the  fuzing  program. 

b.  (S)  The  errors  of  the  fuzing  components  may  be  ascribed  to  the 
error  of  the  integrating  accelerometer  (e)  and  the  error  of  the  timing  device 
(6T).  Provided  that  6T  is  of  the  order  of  20  milliseconds  or  less,  the  timer 
will  have  no  practical  influence  on  the  burst-height  error.  A  value  of  e  of 
0.5  percent  of  the  total  output  of  the  accelerometer  will  give  errors  in  burst 
height  of  the  same  order  of  magnitude  as  the  meteorologically  induced  errors. 

c.  (S)  In  order  to  minimize  the  burst-height  error,  the  mean  atmos¬ 
phere  for  the  location  and  month  should  be  used  in  programming  the  fuzing 
system.  Otherwise  a  bias  may  be  introduced  which  for  some  targets  may  amount 
to  as  much  as  1000  feet  in  altitude  and  one-half  mile  in  range. 

d.  (S)  Based  on  a  precision  for  the  accelerometer  of  0.5  percent  and 
for  the  timer  of  20  milliseconds  or  better,  the  least  dispersion  in  burst  alti¬ 
tude  i3  achieved  when  the  signal  altitude  is  chosen  6000  to  7000  feet  above  the 
burst  altitude.  The  time  of  flight  from  signal  to  burst  with  this  choice  is 
two  to  three  seconds.  This  choice  is  practically  independent  of  season  and 
locale . 

e.  (S)  The  minimum-dispersion  altitude  is  not  sharply  defined.  A 
displacement  of  the  signal  altitude  upwards  from  the  least-dispersion  altitude, 
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by  several  thousand  feet  will  not  increase  the  standard  deviation  of  the  burst- 
height  error  by  more  than  a  few  tens  of  feet. 

f.  (S)  As  a  general  conclusion,  it  can  be  stated  that  a  choice  of 
signal  altitude  anywhere  between  5000  and  15,000  feet  above  the  burst  altitude 
will  give  a  satisfactory  approximation  to  a  least  dispersion. 

g.  (S)  The  adopted  value  of  e  does  influence  these  results.  The 
general  effect  of  an  increase  of  e  is  to  push  the  least-dispersion  signal 
altitude  higher  up. 

h.  (S)  With  a  value  of  £  of  0.5  percent  the  standard  deviation  of  the 
burst-height  errors  will  range  from  about  200  feet  to  700  feet,  depending  on 
season,  location,  and  burst  height. 
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Appendix  A 

(U)  CLOSED-FORM  SOLUTION  FOR  THE  BURST  ERROR 


(U)  The  Re-Entry  Model 

(U)  When  studying  the  atmospheric  perturbations  during  re-entry  it  is  suf¬ 
ficiently  accurate  to  use  an  approximate  re-entry  model.  The  model  used  in  the 
present  study  of  the  re-entry  perturbations  is  given  by  the  equations: 

(A_1)  I?  “  “  ^  sp  w/CjjA  V* 

(A-2)  a  =  ag  =  constant 

(A-3)  cd  “  CD  =  constant 

The  quantities  In  the  equations  above  are : 

V  =  speed  of  vehicle 
t  =  time 

g  =  acceleration  of  gravity 

p  ■  air  density 

W  =  vehicle  weight 

CD  ■  drag  coefficient 

c£  =  hypersonic  drag  coefficient 

A  =  normal  cross  section 

a  =  path  angle 

Og  =  path  angle  at  entry 

(U)  While  equations  (A-l)  through  (A-3)  represent  only  approximate  tra¬ 
jectories,  they  have  sufficient  fidelity  to  be  used  as  basis  for  a  study  of  the 
trajectory  perturbations  caused  by  atmospheric  anomalies. 

(U)  The  method  followed  in  this  study  is  analogous  to  the  method  used  by 
the  author  in  several  other  studies  of  re-entry  problems  [1]  [2]  [3].  The 
object  is  to  obtain  a  closed-form  expression  which  relates  the  re-entry  param¬ 
eter  of  interest  to  the  atmospheric  perturbation  variables.  In  the  case  of  the 
Mark  12  Re-Entry  Vehicle  we  are  interested  in  the  deviations  from  the  pro¬ 
grammed  burst  altitude . 

(U)  The  closed-form  expression  also  points  the  way  to  optimization  of  the 
fuzing  system.  By  proper  selection  of  the  fuzing  parameters  the  burst  disper¬ 
sion  can  be  minimized.  As  it  turns  out  the  selection  is  dictated  by  the 
characteristics  of  the  atmospheric  anomalies  and  by  the  precision  of  the 
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fuzing  components. 

(U)  The  Fuzing  System  In  the  Mark  12  Re-Entry  Vehicle 

(S)  The  Mark  12  fuzing  system  consists  of  an  integrating  accelerometer  and 
a  timing  device.  The  integrating  accelerometer  measures  the  drag  deceleration 
integrated  over  time.  When  a  predetermined  value  of  this  integral  is  reached 
at  an  altitude  we  shall  name  the  signal  altitude  and  denoted  by  subscript  S, 
the  timer  starts  measuring  time,  and  when  a  certain  time  interval  T.  has 
elapsed,  the  burst  mechanism  is  activated.  The  event  of  burst  is  denoted  by 
the  subscript  B. 

(S)  The  purpose  of  this  study  is  to  define  an  optimum  fuzing  system  which 
selects  the  value  of  S  for  a  given  B  such  that  the  error  in  the  burst  altitude 
caused  by  the  combined  effect  of  atmospheric  anomalies  and  component  errors  is 
at  a  minimum. 

(U)  The  Integral  of  Drag  Deceleration 

(S)  Equation  (A-l)  combined  with  the  hydrostatic  equation: 

(A-4)  -gp  =  -  ft 

where  p  is  pressure  and  z  altitude,  and  the  kinematic  relation: 

(A-5)  ||  =  -  V  sin  a 

give 

(A— 6)  V  =  VE  exp  (-p/2P) 

VE  is  speed  at  entry  of  the  atmosphere  where  p  **  0.  The  constant 
(A-7)  P  =  sin  a  W/CjjA 

which  has  the  dimension  of  pressure,  represents  the  atmospheric  pressure  at 
peak  deceleration.  For  the  Mark  12  Re-Entry  Vehicle  W/CpA  =  2000  lbs  per 
square  foot,  and  with  a  =  20*,  2P  becomes  655  mb. 

(U)  The  output  of  the  integrating  accelerometer  from  a  threshold  value, 
which  for  practical  purposes  may  be  chosen  as  zero,  to  the  signal  altitude  S, 
may  now  be  written,  observing  (A-6): 

tS 

(A-8)  I  -  f  fjf  dt  =  VE  -  Vs  [1  -  exp  (-Pg/2P )  ] 
tE 
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(U)  It  will  be  seen  from  the  preceeding  equations  that  the  gravity  accel¬ 
eration  g  sin  a  has  been  neglected  beside  the  drag  decelerations.  However,  as 
we  are  interested  in  the  environmental  perturbations  this  is  of  little  conse¬ 
quence,  as  the  gravity  field  remains  unaffected  by  the  perturbations. 

(U)  Time  of  Flight  from  Signal  to  Burst 

(U)  Returning  to  equations  (A-l)  and  (A-6)  we  may  write  for  the  time  of 
flight  T  between  the  output  signal  at  S  to  the  burst  at  B: 

B  PB  PB 

(A-9)  T  =/dt  -  $  -g  VE  Sin  a/  0  exp  <p/2P)  dln  P 

S  Ps  Ps 

(U)  In  (A-9)  use  has  been  made  of  the  equation  of  state  for  moist  air: 

(A- 10)  P  = 

where  9  Is  the  virtual  temperature  and  R  the  gas  constant  for  dry  air. 

(U)  We  now  Introduce  a  new  non-dimensional  pressure  variable,  cd,  defined 
by: 

(A-ll)  dou  =  exp  (p/2P)  din  p 


or,  integrated: 

p/2P 

co  =  J'  exp  §  din  §  =  Ei  (p/2P) 


where  Ei(x)  =  J'  exp  §  din  ?  is  the  exponential  integral. 


(U)  Equation  (A-9)  may  now  be  written: 


“b 

T  '  g  V,  sin's/  9 

“s 

The  limits  of  the  Integral  are  determined  from  cDg 


E 1  (Pq/2P)j  cdb  =E i  (Pb/2P). 


(U)  The  relationship  in  the  atmosphere  between  temperature  9  and  pressure 
p  or  cd  is  observed  by  means  of  radiosondes  and  rocketsondes .  The  functional 
relationship  is  usually  recorded  in  the  form  of  a  trace  in  some  aerologlcal 
diagram,  and  it  has  been  found  convenient  to  develop  a  graphical  means  of 
evaluating  this  trace  in  terms  of  the  re-entry  performance  of  vehicles.  The 
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so-called  Theta-Omega  Diagram  described  in  [1]  has  been  constructed  for  this 
purpose. 

(U)  The  Perturbation  Equation 

(u)  A  re-entry  system  and  its  respective  fuzing  system  is  programmed  to 
perform  in  a  certain  standard  atmosphere .  For  our  purpose  it  is  not  necessary 
at  this  junction  to  prescribe  which  standard  is  to  be  used.  It  may,  for  in¬ 
stance,  be  the  1962  U.S.  Standard  Atmosphere  [4],  or  some  climatological 
profile  valid  for  a  specific  month  and  location. 

(U)  The  fuzing  program  is  then  based  on  zero  errors  of  the  accelerometer 
and  the  timing  device.  Quantities  referring  to  this  idealized  behavior  in  a 
standard  atmosphere  and  with  no  component  error  will  be  denoted  by  an  asterisk. 

(U)  In  the  real  case  the  integrating  accelerometer  will  have  an  error  61 
defined  by: 

(A-13)  I  =  I*  +  61 

where  I*  is  the  programmed  value.  In  the  following  it  Is  assumed  that  61  is 
roughly  proportional  to  I*: 

(A-14)  61  -  el* 

so  that  e  is  a  constant  which  characterizes  the  instrument. 

(U)  The  error  of  the  timing  device  6T  is  defined  by: 

(A-15)  T  =  T*  +  6T 

where  T*  is  the  programmed  time. 

(U)  The  accelerometer  error  61  causes  a  deviation  in  the  pressure  at  sig¬ 
nal  altitude  which  may  be  computed  from  equation  (A-8).  By  a  variation  of 
equation  (A-8)  and  using  (A-14),  we  find: 

(A-l6)  6pg  =  2P  [exp  (p|/2P)  -  1]  e 

(U)  In  equation  (A-16)  p|  is  the  pressure  at  the  programmed  signal  alti¬ 
tude,  ps,  the  corresponding  value  in  the  real  case,  and  6pg  is  defined  by: 

(A-17)  PS  -  Pg  +  ^Pg 

(U)  Equations  (A-15)  and  (A-12)  combine  to  give: 
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B  g  V-.  sin  a 

9*  do) - ^ -  6T  =  0 

Pg  is  the  programmed  pressure  at  burst,  9*  the  virtual  temperature  in  the  stand¬ 
ard  atmosphere,  while  the  quantities  without  the  asterisk  again  denote  the  real 
case . 

(U)  We  further  define: 

(A-19)  0  =  9*  +  69 

and 


PB  P 

(A-18)  f  6  to-/ 
PS  PS 


(A-20)  Pg  =  Pg  +  6pB 


(u)  Expanding  (A-l8)  and  using  (A- 19)  and 
Pg 

g  VE  sin 

60  do>  +  0*  6a-  h 


(A— 2! )  £ 


SB  0<UB  "  ®S  IT 


PS 


(A-20)  we  obtain: 
a 

-  6T  =  0 


where  we  have  retained  first  order  terms  only,  making  use  of  such  high-accuracy 
approximations  as  the  following: 


®B  «®B 

6S  6“S 


(U)  6tOg  and  6<0g  are  computed  from  (A-ll): 

6o>s  =  exp  (p£/2P) 

ps 

and  from  (A-l6): 

(A-22)  6cds  =  2P  exp  (p^/2P)  [exp  (p|/2P)  -1]  e/pg 
(A-23)  6u)g  =  exp  (p£/2P)  ipg/Pg 

(U)  The  pressure  deviation  at  burst,  6pB,  is  related  to  the  altitude 
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deviation  6zg  through  the  hydrostatic  equation  (A-4).  We  obtain: 

R  Q* 

(A- 24)  6zb - 1  6pB  +  Db 

g  PB 

(U)  In  the  last  equation,  Dg  is  the  so-called  D-value  at  burst  altitude. 
The  D-value  is  defined  as  the  altitude  of  a  given  pressure  in  the  real  atmos¬ 
phere  minus  the  altitude  of  the  same  pressure  in  the  standard  atmosphere.  The 
D-value  appears  on  the  right-hand  side  of  equation  (A-24)  because  the  geometric 
measurement  of  z  refers  to  the  ground  as  a  reference,  whereas  pressure  p  does 
not  make  use  of  such  a  reference. 

(U)  (A-24)  together  with  (A-22)  and  (A-23)  inserted  into  (A-21)  gives  us 
our  final  equation: 

PB 

(A-25)  6zB  =  Dg  +  |  exp  (-p^/2P)  J  <50  doj 


R  9  c 


g  Pc 


exp 


-  Vg  sin  a  exp 


6T 


1 


2Pe 


(U)  (A-25)  relates  the  deviation  in  the  burst  altitude,  6zg,  to  two  mete¬ 
orological  variables  plus  two  fuzing  component  errors. 

a.  (U)  The  meteorological  variables  are: 

(l)  (U)  The  mean  temperature  deviation  between  signal  and  burst. 


represented  by  the 


(2)  (U)  The  D-value  at  burst,  Dg. 

b.  (U)  The  fuzing  component  errors  are: 

(1)  (U)  The  relative  error,  e,  of  the  output  from  the  integrating 
accelerometer. 

(2)  (U)  The  timing  error  6T  of  the  device  which  counts  time  between 
signal  and  burst.  [Errors  of  the  detonation-activation  mechanism,  etc.,  have 
been  lumped  with  the  timing  error.] 
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(U)  Once  the  signal  altitude  and  the  burst  altitude  have  been  determined 
the  somewhat  involved  multipliers  for  the  variables  become  constants,  and  we 
can  write: 


(A-26)  6zB  «  Dg  +  M  -  Ne  -  Q  6T 
where 


(A-27)  66 


—‘I  > 

-  to* 


66  dto 


“B  -  “S  pi 

equals  the  mean  temperature  deviation  from  signal  to  burst, 


M  =  -  exp  (-pg/2P)  {tog  -  tog) 


2p  r  e! 


g  Pc 


(-^)h(§)-1] 

(-3) 


exp 


Q  =  V„  sin  a  exp 


Q  depends  on  burst  altitude  only,  M  and  N  on  both  signal  altitude  and  burst 
altitude . 

(U)  In  the  following  the  1962  U.S.  Standard  Atmosphere  [4]  will  be  used  as 
a  reference  atmosphere.  The  multipliers  M,  N,  and  Q  are  based  on  this  atmos¬ 
phere  and  the  D-values  also  refer  to  this  standard  atmosphere. 

(U)  Figures  7a  through  8b  show  the  mean  and  the  standard  deviation  of  the 
D- value  at  sea  level  for  the  months  of  January  and  July.  Figures  9  and  10  show 
the  corresponding  values  at  18,000  feet.  For  burst  altitudes  between  sea  level 
and  18,000  feet  the  mean  and  the  standard  deviation  of  the  D-value  may  be  ob¬ 
tained  by  a  linear  Interpolation. 

(U)  M  has  been  graphed  in  Figure  4.  The  units  are  in  feet  per  degree 
Celsius. 

(U)  N  is  graphed  in  Figure  5.  Units  are  feet  per  percent  relative  accel¬ 
erometer  error. 

(S)  Q  is  depicted  in  Figure  6  in  feet  per  millisecond  of  error  of  the 
timer.  As  6T  is  of  the  order  of  one  millisecond  in  the  Mark  12  Re-Entry  Vehicle, 
it  is  immediately  obvious  that  the  timing  error  contributes  only  insignificantly 
to  the  over-all  error. 

(S)  If  we  chose  e  =  .5  x  10~8  as  a  representative  root-mean- square  value 
of  the  accelerometer  error,  it  is  seen  that  the  D-value  and  the  contributions 
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of  the  accelerometer  error  and  the  temperature  deviation  are  of  comparable  mag¬ 
nitude  . 

(U)  Of  the  various  error  sources  only  the  meteorological  variables  are 
correlated,  as  there  seems  to  be  no  reason  to  assume  that  the  accelerometer 
error  or  the  timing  error  depend  on  each  other  or  on  the  environment.  These 
relations  are  explored  further  in  Appendix  B. 
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Appendix  B 

(U)  STATISTICAL  DISTRIBUTION  OF  THE  BURST- HEIGHT  ERROR 

(U)  Equation  (A-26)  of  Appendix  A  gives  the  burst-height  error  dZg  as  a 
linear  combination  of  four  error  sources. 

(A-26)  dzB  -  DB  +  M  6§  -  Ne  -  Q  6T 

(U)  For  a  given  burst  altitude  and  a  given  signal  altitude  M,  N,  and  Q  are 
constants. 

(U)  Dg  is  the  D-value  at  the  nominal  burst  altitude. 

(U)  68  is  the  mean  temperature  deviation  from  the  nominal  value  in  the 
layer  from  signal  to  burst. 

(U)  e  is  the  relative  error  of  the  integrating  accelerometer.  In  the  fol¬ 
lowing  calculations  a  rms-value  of  e  =0.5  x  10"*2  will  be  used. 

(U)  6T  is  the  error  of  the  timing  device.  A  rm3-value  of  6T  «  10-a  second 
will  be  used  in  the  following. 

(U)  We  will  assume  normal  distributions  of  Dg,  68,  e,  and  6T,  with  zero 
means  for  €  and  6T 

(U)  If  the  monthly  climatological  average  atmosphere  has  been  used  in  the 

A 

programming  of  the  fuzing  system,  the  means  of  Dg  and  60  are  also  zero,  and  we 
are  only  concerned  with  the  dispersion  of  these  quantities  about  their  clima¬ 
tological  averages. 

;)  Of  these  four  error  sources  only  the  atmospheric  parameters  are  likely 
to  have  significant  correlations. 

(U)  In  the  following,  a(x)  denotes  the  standard  deviation  of  x: 

(B-l)  ff(x)  =  [(x  -  x)8]^ 

where  the  bar  indicates  averaging  over  the  sample,  and  r(x,y)  denotes  the 
linear  correlation  between  x  and  y. 

(U)  With  these  conventions  we  obtain  from  (A-26)  for  the  variance  of  6z0: 
(B-2)  a2(6zg)  =  cr*(Dg)  +  M*a*(0)  +  N*a*(e)  +  Q*aa(6T)  +  2Md(Dg)  a(8)  r(Dg,0) 

(U)  The  D-Value  at  Burst  as  an  Error  Source 

(U)  Figures  7,  8,  9,  and  10  show  the  global  distribution  of  ff(Dg)  at  sea 


15 


Technical  Report  178 


September  1964 


level  and  at  18,000  feet  for  a  winter  month  (January)  and  a  summer  month  (July) 
A  linear  interpolation  with  altitude  will  give  a  good  approximation  to  the 
value  at  a  burst  altitude  below  18,000  feet. 

(U)  Figure  11  shows  o(Dg)  as  a  function  of  “Ititude  for  two  stations  in 
the  U.S.S.R.  during  summer  and  winter. 


(U)  Anomaly  of  Mean-Layer  Temperature  from  Signal  to  Burst  as  an  Error  Source 

(U)  Meteorological  archives  do  not  as  a  rule  contain  statistics  on  mean- 
layer  temperature,  but  contain  statistics  on  point  values  and  correlations 
between  point  values.  This  still  permits  us  to  calculate  the  wanted  distribu¬ 
tions.  The  methods  used  are  indicated  in  the  following. 

(U)  From  equation  (A-27)  it  follows  that: 


O). 


(B-3)  aa(e)  =  (cug  -  a >s)~a  f  du  f  o(S?)  r(0^,0?)  dC 

“s  “s 

where  .ie  variables  p.  and  £  are  dummies  for  o>. 


ox 


B 


(U)  For  the  finite-difference  evaluation  of  (B-3)  we  proceed  as  follows: 
Let  the  Interval  from  signal  to  burst  contain  n  data  points,  of  which  the  sig¬ 
nal  and  burst  altitudes  are  two  points,  as  indicated  in  Figure  12.  To  each 
point  is  ascribed  a  value  of  co  and  a(0);  and  the  correlation  coefficient 
r(01,6j)  for  any  pair  of  points  is  known. 

(U)  When  each  temperature  sounding  in  the  sample  is  approximated  by 
straight-line  segments  between  the  tabulated  data  points,  we  obtain  for  the 
distribution  (B-3): 

n  n 

(B-4)  aa(0)  =  (cog  -  ojg)"*  ^  ^  o(61)  o(9j)  r(91,ej)  k±kj 

1=1  j=l 

where  the  K's  are  given  by: 


kx  =  (o>2  -  m^/2 
kn  =  K  -  “n-l^2 

k±  =  (<o1+1  -  <o1_1)/2;  i  =  2,  3,  . . .  (n-1) 

(U)  Figure  13  shows  M8aB(0)  for  the  two  stations  in  U.S.S.R.  during  Janu¬ 
ary  and  July.  Figure  11  shows  a(0)  as  a  function  of  altitude. 


(U)  The  Integral  Accelerometer  as  an  Error  Source 
(U)  We  will  assume 
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ff(«)  »  ,5x  10“ 2 
and  further 


r(DB, e)  =  r(@, c)  =  r(6T,e)  =  0 

(U)  Figure  5  shows  twice  the  error  due  to  this  effect. 

(U)  The  Timing  Device  as  an  Error  Source 
(U)  A  value 

a(6T)  =  10“ 3  second 

gives  the  contribution  illustrated  in  Figure  6. 

(U)  Errors  Due  to  Correlation  of  Mean  Temperature  and  D-Value 

(U)  The  last  term  in  equation  (B-2)  contains  r  (Dg, 9)  which  has  to  be 
evaluated  by  special  technique  due  to  the  same  reasons  as  indicated  for  c(0). 

(U)  From  the  definition  of  the  vertical  average,  given  In  Appendix  A: 


“b 

x  =  (cuB  -  (Ug)"1  J  x 


dtu 


combined  with  the  classical  definition  of  the  correlation  coefficient 
r(x,y)  =  [a(x)  a(y)]"1  xy 
we  obtain,  after  some  calculation: 


B 


(U)  Using  the  same  finite-difference  scheme  as  before  (Figure  12^  we  can 


(B-5) 

p(Db,0)  = 

(u) 

Using  the 

write: 

(B— 6 ) 

r(Dg,9)  = 

(u) 

Figure  14 

dec 


1=1 


locations  in  summer  and  winter. 
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(U)  The  Dispersion  of  the  Burst  Height 

(U)  Prom  each  of  the  error  sources  discussed  above,  we  may  construct  the 
distribution  of  the  total  error  and  Its  dependence  on  the  altitudes  of  signal 
and  burst.  Figures  1  and  2  show  the  result  of  our  calculations  for  the  four 
atmospheres  chosen  as  examples .  In  the  six  panels  of  Figure  1,  panels  A,  B, 

C,  D,  and  E  show  the  contributions  from  each  of  the  terms  in  equation  (B-2)  to 
the  variance  of  6zg.  Panel  F  shows  the  standard  deviation  of  the  burst  error 
and  represents  the  final  result  of  our  calculations. 

(U)  Figure  2  shows  the  end  result  for  all  four  atmospheres.  Table  1  con¬ 
tains  the  results  in  tabular  form. 
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Contribution  to  Variance  [Columns  A,  B,  C,  and  D  x  10-4  (feet)2]  and  Standard  Deviation  of  the  Burst-Height  Error 
[F  (feet)],  (a)  Station  38457,  July;  (b)  38457,  January;  (c)  Station  23552,  July;  (d)  23552,  January;  S  is  signal 
altitude;  B  is  burst  altitude. 
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STANDARD  DEVIATION  OF  BURST  ERROR,  FEET 
(FOUR  YEARS  DATA:  1955-1958) 
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13.4  11.82  10.38  9.06  7.71  6.42 


BURST  ALTITUDE,  THOUSANDS  OF  FEET 
LAG, SIGNAL  TO  BURST,  IN  SECONDS 

MINUMUM  DISPERSION  AGAINST  TIME  OF  FLIGHT 

Figure  3.  (S) 
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Figure  5.  (S) 
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Figure  6.  (3) 
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Figure  7a.  (U) 


26 

SECRET 

(This  page  Is  unclasslfl< 


September  1964 


Technical  Report  178 


Technical  Report  178 


September  1964 


V 


Technical  Report  178 


September  1964 


Figure  9.  (U) 
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Figure  10.  (U) 
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Figure  12.  (U) 
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DEPARTMENT  OF  THE  AIR  FORCE 

Air  Force  Combat  Climatology  Center  (AFWA) 
Asheville,  north  Carolina  28801 -5002 


4  April  2006 


MEMORANDUM  FOR  DTIC-OQ 

ATTENTION:  LARRY  DOWNING 

8725  JOHN  J.  KINGMAN  ROAD 
FORT  BELVOIR,  VA  22060-62 1 8 

FROM:  Air  Force  Weather  Technical  Library 
151  Patton  Ave,  Rm  120 
Asheville,  NC  28801-5002 

SUBJECT:  CHANGE  CLASSIFICATION  AND  DISTRIBUTION  STATEMENTS 


1 .  AD353247  -  Least  dispersion  of  the  Mark  1 2  re-entry  vehicle,  Sep  1 964. 

2.  AD8 19674  -  Climatic  data  summaries  for  selected  Asian  and  Pacific  stations;  Republic  of  Vietnam  stations,  Apr 
1967. 

3.  AD256922  -  Bibliography  of  documents  prepared  by  weather  staff  sections  of  Headquarters  Army  Air  Corps  and 
Headquarters  Army  Air  Forces,  1937-June  1945,  Apr  1961. 

4.  AD2597 1 0  -  Technical  publications  and  documents  of  Air  Weather  Service  field  activities,  Part  one:  1 945-6 1  with  a 
chronology  of  AWS  organizations,  1945-61,  June  1961. 

5.  AD156453  -  Charts  of  maximum  and  minimum  thickness  lines  1000/500  MB,  Northern  hemisphere,  1946-1950, 
Sep  1952. 

6.  AD244550  -  Technical  publications  and  documents  of  Air  Corps  and  AAF  Field  weather  activities,  Part  two:  1937- 
1945  with  a  chronology  of  AC  and  AAF  Weather  regional  organizations,  1937-1945,  Sep  1960. 

All  the  above  documents  need  to  be  changed  to  "Approved  for  Public  Release,  Distribution  Unlimited”  please. 


SUSAN  A.  TARBELL 
Librarian,  Classified  Custodian, 
DTIC  Point  of  Contact 
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